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We  present  a  dissipative  particle  dynamics  simulation  study  on  the  formation  of  nanostructures  of  symmetric  diblock 
copolymers  confined  between  planar  surfaces  with  and  without  nanopatterns.  The  nanopatterned  surface  is  mimicked 
by  alternating  portions  of  the  surface  that  interact  differently  with  the  diblock  copolymers.  The  formation  of  the 
diblock— copolymer  nanostructures  confined  between  the  planar  surfaces  is  investigated  and  characterized  by  varying 
the  separation  width  and  the  strength  of  the  interaction  between  the  surfaces  and  the  diblock  copolymers.  For  surfaces 
with  nanopatterns,  we  also  vary  both  the  mutual  area  and  location  of  the  nanopatterns,  where  we  consider  nanopatterns 
on  the  opposing  surfaces  that  are  vertically  (a)  aligned,  (b)  staggered,  and  (c)  partially  staggered.  In  the  case  of  planar 
slits  without  nanopatterns,  we  observe  the  formation  of  perpendicular  and  parallel  lamellar  phases  with  different 
numbers  of  lamellae.  In  addition,  the  symmetric  diblock  copolymers  self-assemble  into  adsorbed  layer  and  adsorbed 
layer— parallel  lamellar  phases  and  a  mixed  lamellar  phase  when  the  opposing  surfaces  of  the  planar  slits  are  modeled  by 
different  types  of  wall  beads.  In  the  case  of  nanopatterned  planar  slits,  we  observe  novel  nanostructures  and  attempt  to 
rationalize  the  diblock  copolymer  self-assembly  on  the  basis  of  the  behavior  that  we  observed  in  the  planar  slits  without 
nanopatterns.  In  particular,  we  investigate  the  applicability  of  predicting  the  structures  formed  in  the  nanopatterned  slits 
by  a  superposition  of  the  observed  structures  in  slits  without  nanopatterns. 


1.  Introduction 

Diblock  copolymers  consist  of  two  chemically  incompatible 
blocks  A  and  B  that  are  connected  end-to-end.  When  diblock 
copolymers  are  quenched  under  the  order— disorder  temperature, 
the  A  and  B  blocks  tend  to  form  separated  phases;  however,  they 
are  prevented  from  doing  so  by  the  end-to-end  connectivity. 
Instead,  the  diblock  copolymers  self-assemble  into  nanostructures 
that  minimize  contacts  between  the  A  and  B  blocks.  The  mor¬ 
phology  of  the  nanostructures  depends  on  the  length  ratio  of  the 
two  blocks  and  the  strength  of  chemical  incompatibility  between 
the  blocks.  The  natural  tendency  of  diblock  copolymers  to  self- 
organize  on  the  mesoscale  into  dense,  regular  arrays  of  nano¬ 
structures  makes  them  attractive  for  various  nanotechnology- 
related  applications.1  In  the  bulk  phase,  diblock  copolymer 
systems  form  four  classical  microphases:  the  lamellar  phase  for 
symmetric  (where  the  lengths  of  the  A-  and  B-blocks  are  equal)  or 
slightly  unsymmetric  diblock  copolymers  and  the  cylindrical, 
gyroid,  and  micellar  phases  otherwise.2  However,  in  confinement, 
adsorption  interactions,  symmetry  breaking,  structural  frustra¬ 
tion,  confinement-induced  entropy  loss,  and  heterogeneous  sur¬ 
faces  play  a  dominant  role  leading  to  nanostructures  that  differ 
front  those  found  in  the  bulk.  From  a  technological  viewpoint,  the 
intriguing  prospects  of  diblock  copolymer  self-assembly  in  nano¬ 
patterned  confinement  are  indeed  the  novel  nanostructures  that 
can  form,  but  the  possibility  that  they  may  serve  as  templates  for 
other  types  of  nanostructures  is  also  intriguing. 
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Our  understanding  of  the  self-assembly  of  symmetric  diblock 
copolymers  confined  between  planar  slits  without  nanopatterns 
continues  to  deepen  via  experimental,  theoretical,  and  simulation 
studies;3,4  for  example,  the  influence  of  copolymer  polydispersity 
was  recently  studied.5  For  these  systems,  typically  each  slit  wall 
has  a  substantial  preferential  affinity  toward  one  of  the  blocks, 
which  results  in  lamellae  aligning  in  a  parallel  orientation  with 
respect  to  the  walls  as  shown  in  Figure  la.  Ideally,  the  parallel 
lamellar  phase  adopts  lamellae  of  thickness  vUlam,  where  v  is  an 
integer  corresponding  to  the  number  of  lamellae  and  W®AM  is  the 
bulk-phase  lamellar  spacing.  In  cases  where  the  natural  spacing  of 
the  lamellae  is  incommensurate  with  the  separation  of  the  walls  or 
slit  width  W(i.e. ,  Wfv  W°AM  for  any  value  of  v),  a  parallel  phase 
would  be  too  energetically  and  structurally  frustrated.  Thus,  to 
avoid  this  situation  the  lamellae  prefer  to  align  perpendicular  to 
the  walls,  as  shown  in  Figure  lb.  Whereas  a  perpendicular 
lamellar  phase  cannot  accommodate  the  energetically  favorable 
layering  of  the  preferentially  adsorbed  block,  its  lamellae  spacing 
is  free  to  relax  to  its  bulk-phase  spacing.  Lamellae  can  also 
orientate  to  a  mixed  lamellar  phase,  as  shown  in  Figure  lc,  with 
v  parallel  lamellae  adsorbed  at  one  surface  and  perpendicularly 
orientated  lamellae  at  the  other  surface.  Although  the  presence  of 
an  interface  between  the  two  lamellae  orientations  is  energetically 
costly,  this  cost  is  compensated  for  by  the  adsorbed  layer  at  the 
wall.  A  mixed  lamellar  phase  is  slightly  unstable  with  respect  to  a 
perpendicular  lamellar  phase;  however,  this  morphology  can  be 


(3)  Binder,  K.  Adv.  Polym.  Sci.  1999,  138,  1-89. 

(4)  Fasolka,  M.  J.;  Mayes,  A.  M.  Annu.  Rev.  Mater.  Res.  2001,  31,  323-355. 

(5)  Han,  Y.;  Cui,  J.;  Jiang,  W.  Macromolecules  2001,  41,  6239-6245. 


Langmuir  XXXX,  XXX (XX),  XXX-XXX 


DOI:  10.1021/la903200j  A 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

03  OCT  2009  2- REPORT  TYPE 

3.  DATES  COVERED 

00-00-2009  to  00-00-2009 

4.  TITLE  AND  SUBTITLE 

Self-Assembly  of  Symmetric  Diblock  Copolymers  in  Planar  Slits  with  and 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

E.  Hala  Laboratory  of  Thermodynamics, Institute  of  Chemical  Process 
Findamentals  of  ASCR,v.  v.  i.  Rozvojova  135,165  02  Prague  6-Suchdol, 
Czech  Republic, 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

18.  NUMBER  19a.  NAME  OF 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Public  Release 

unclassified  unclassified  unclassified 

15 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Downloaded  by  Martin  Lisal  on  October  19,  2009  I  http://pubs.acs.org 
Publication  Date  (Web):  October  19,  2009  I  doi:  10.1021/la903200j 


Article 


Petrus  et  al. 


Figure  1.  Examples  of  simulation  configurations  for  (a)  parallel,  (b)  perpendicular,  and  (c)  mixed  lamellar  phases  obtained  by  dissipative  particle 
dynamics  simulations  of  symmetric  diblock  copolymers  A5B5  in  planar  slits  without  nanopattems.  Green  and  gold  spheres  represent  A  and  B  beads, 
respectively,  whereas  gray  spheres  represent  wall  beads.  B  beads  have  preferentially  affinity  for  the  walls  and  thus  are  adsorbed  onto  the  walls. 


preserved  by  imposing  a  small  asymmetry  between  the  copolymer 
block  lengths.6’7 

Similarly,  many  experimental,  theoretical,  and  simulation 
studies  have  been  performed  on  the  self-assembly  and  adsorption 
of  copolymers  on  nanopatterned  surfaces;  however,  we  are  not 
aware  of  a  systematic  study  involving  symmetric  diblock  copoly¬ 
mers.  A  few  notable  computational  studies  are  the  work  of  Chen 
et  al.,8’9  where  Monte  Carlo  simulations  of  the  self-assembly  of 
multiblock  copolymers  on  nanopatterned  surfaces  were  carried 
out.  Also  worth  noting  are  the  works  of  Genzer10  and  Semler  and 
Genzer,* 11  where  a  three-dimensional  self-consistent-field  model 
was  used  to  study  the  adsorption  of  copolymers  on  planar 
substrates  composed  of  two  chemically  distinct  regions.  These 
studies,  as  well  as  the  recent  experimental  nanoassembly  of 
symmetric  diblock  copolymers  on  chemically  striped  patterned 
surfaces,12  demonstrate  the  strong  influence  of  nanopatterns  on 
the  copolymer  conformational  and  adsorption  behavior  as  well  as 
the  ability  of  the  copolymers  to  recognize  the  surface  motif,  which 
also  contributes  to  the  self-assembly. 

In  this  work,  we  use  the  dissipative  particle  dynamics  (DPD) 
method  to  explore  the  spontaneous  nanostructure  formation  of 
symmetric  diblock  copolymers  confined  between  planar  surfaces 
both  with  and  without  nanopattems.  The  surfaces  are  modeled  as 
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dense  layers  of  DPD  beads,  where  the  nanopattems  are  modeled 
by  specifying  alternating  portions  of  the  surface  to  interact 
differently  with  the  diblock  copolymers.  The  competition  between 
the  parallel  and  perpendicular  lamellar  alignments  that  occurs  as  a 
consequence  of  the  nanopattems  leads  to  novel  nanostructures. 
The  self-assembly  of  nanostructures  in  planar  slits  is  investigated 
and  characterized  by  varying  the  slit  width  and  the  interaction 
between  the  walls  and  the  diblock  copolymers.  For  nanopatterned 
surfaces,  we  also  varied  both  the  mutual  area  and  location  of  the 
nanopattems,  where  we  considered  nanopattems  on  the  opposing 
surfaces  that  are  vertically  (a)  aligned,  (b)  staggered,  and  (c) 
partially  staggered.  Qualitative  phase  diagrams  are  presented  for 
the  various  systems  along  with  the  corresponding  conformational 
behavior  of  the  diblock  copolymers.  We  make  no  attempt  to  locate 
the  precise  microphase  separation  boundaries;  rather  we  attempt 
to  quantify  the  various  effects  on  self-assembly. 

The  remaining  sections  of  the  article  are  organized  as  follows. 
The  simulation  methodology  is  described  in  section  2,  along  with 
the  DPD  model  of  the  confined  diblock  copolymer  system  and 
computational  details.  Results  are  discussed  in  section  3.  Finally, 
section  4  gives  our  conclusions.  Note  that  throughout  the  article 
we  use  the  terms  surface  and  wall  interchangeably  when  describing 
the  behavior  in  the  planar  slits. 

2.  Simulation  Methodology 

2.1.  DPD  for  Confined  Diblock  Copolymer  Systems. 

In  our  DPD  simulations,13’14  the  diblock  copolymer  chain  is 
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modeled  as  a  collection  of  point  beads  that  represent  lumps  of  the 
chain  containing  several  segments.15  Confining  walls  are  modeled 
by  creating  dense  layers  of  DPD  beads  that  represent  several 
atoms,  where  the  DPD  beads  can  be  arranged  either  randomly16 
or  on  a  lattice. 1 7  The  density  of  the  wall  layers  must  be  sufficiently 
large  to  avoid  penetration  of  the  diblock  copolymer  beads  into 
these  confining  walls. 

In  this  study,  all  DPD  beads  are  defined  by  a  mass  m„  position 
r„  and  velocity  \h  and  interact  with  each  other  via  a  pairwise,  two- 
body,  short-ranged  force  F  that  is  written  as  the  sum  of  a 
conservative  force  Fc,  dissipative  force  F°,  and  random  force  FR: 


+  F"  +F* 


(1) 


For  the  linear  diblock  copolymer  chains,  Fc  includes  a  soft 
repulsion  force,  F  ,  acting  between  two  nonbonded  beads  and  a 
harmonic  spring  force,  FCs,  acting  between  bonded  beads  in  a 
chain.  The  DPD  diblock  copolymer  chains  are  flexible  because  no 
additional  constraints  such  as  bond  bending  or  bond  torsion  are 
included.  For  a  wall  bead,  Fc  contains,  in  addition  to  FCr,  a  wall 
force  FCw  that  tethers  a  wall  bead  to  its  initial  position  rw. 
Expressions  for  FCr,  FCs,  and  FCw  are  given  by 


F 


c 

ij 


t  ij  <  rc 

rij—rc 


(2) 


F<=s  =  -Ks(nj-r0)rf 
rij 

and. 

Ffw  =  — £w(r,—  rw>/) 


(3) 

(4) 


respectively.  In  eqs  2—4,  ay  is  the  maximum  repulsion  between 
bead  i  and  bead  j,  ty  =  r,-  —  r,,  ty  =*•  jr^l,  rc  is  the  cutoff  radius,  Ks  is 
the  spring  constant,  r0  is  the  equilibrium  spring  length,  and  Kw  is 
the  wall-bead  tethering  constant. 

The  remaining  two  forces  F°  and  FR  are  given  by 

f  (5) 

(6) 

where  aP(ry )  and  a)R(ry)  are  the  weight  functions  that  vanish  for 
Ty  >  rc,  Yy  is  the  friction  coefficient,  Oy  is  the  noise  amplitude,  Sy  = 
v,;  -  v,-,  is  a  Gaussian  random  number  with  a  zero  mean  and  a 
unit  variance  that  is  chosen  independently  for  each  pair  of 
interacting  beads,  and  A?  is  the  time  step. 

Espanol  and  Warren18  showed  that  the  system  samples  the 
canonical  ensemble  and  obeys  the  fluctuation— dissipation 
theorem  (in  the  limit  of  At  — *  0)  if  the  following  relations  hold 
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where  T  is  the  temperature  and  kB  is  the  Boltzmann  constant. 
aiD(r,y)  and  ajR(ry )  are  typically  chosen15  as 


«D(r;;)  =  [ft»R(n/)]“ 


rij  <  rc 

rij—rc 


(9) 


The  evolution  of  DPD  beads  in  time  t  is  governed  by  Newton’s 
equations  of  motion 


dr; 
d  t 


=  v;« 


(10) 


«,^=f;(0=E(F'+F°+F^) 

¥J 

Our  mesoscale  model  of  diblock  copolymer  A„B„,  consists  of 
A„  and  B,„  blocks  connected  end  to  end  where  n  and  in  are  the 
number  of  A  and  B  beads,  respectively.  To  determine  the  values  of 
ay,  we  used  expressions  obtained  by  Groot  and  Warren  that  map 
DPD  polymer  models  onto  Flory— Huggins- type  models.15 
Groot  and  Warren  derived  the  following  expressions  for  the  like 
interactions 


flAA^c  _  flBB'c  _  75 
A'b T  ~  kBT  p/'c3 

and  for  the  unlike  interactions 

flAB rc  _  «BBfi:  ,  «2 

kBT  kBT  a\  (, n+m)a 3 

where  p  is  the  bead  number  density,  ^ab  is  the  effective  Flory— 
Huggins  interaction  parameter,  and  au  a2,  and  a3  are  constants 
dependent  on  p/-c3  (e.g.,  ci\  =  3.27,  a2  =  3.9,  and  a3  =  0.51  when 
p/-c3  =  3, 19  which  is  used  in  this  work). 

The  confining  walls  are  composed  of  four  layers  of  DPD  beads 
arranged  on  a  face-centered-cubic  lattice  (following  the  work  of 
Malfreyt  and  Tildesley17)  with  a  wall  density  of  approximately 
3.5p.  Although  not  presented  here,  the  calculated  A-  and  B-bead 
density  profiles  confirmed  that  the  wall  density  used  avoids  the 
penetration  of  the  diblock  copolymer  beads  into  the  confining 
walls.  For  slits  without  nanopattems,  the  opposing  walls  are 
either  created  from  DPD  beads  of  the  same  type,  nq,  or  created 
from  DPD  beads  of  different  types,  w\  and  w2,  which  possess 
preferential  affinity  toward  the  B  and  A  beads,  respectively.  For 
slits  with  nanopatterns,  each  wall  contains  alternating  areas  or 
strips  created  from  nq  and  w2  beads.  Because  the  wall  beads  are 
tethered  to  their  lattice  positions  by  stiff  harmonic  springs,  no 
additional  repulsion  interactions  between  the  wall  beads  them¬ 
selves  are  needed;  therefore,  we  set  aWiWi  =  an.2W ,  =  aWiWl  =  0. 
Furthermore,  for  the  sake  of  simplicity,  the  repulsion  interactions 
between  A  and  w2  beads  and  between  B  and  vtq  beads  were 
considered  to  be  the  same  as  the  like  interactions  for  diblock 
copolymers  (i.e.,  we  set  aAwi  =  aBn.t  =  aAA  =  Abb)-  Values  of  aAW[ 
and  «b„2  were  then  varied  in  this  study. 

2.2.  Computational  Details.  Throughout  this  article,  we 
used  the  following  reduced  units;  rc  is  the  unit  of  length,  kBT  is 
the  unit  of  energy,  and  the  mass  of  a  DPD  bead  is  the  unit  of  mass. 
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Using  these  reduced  units,  we  set  p  =  3,  07  =  3,  Kn.  =  100  and 
At  =  0.03  and  assumed  that  the  DPD  beads  corresponding  to 
the  diblock  copolymers  and  the  walls  have  the  same  mass, 
/nDPD.  We  investigated  the  behavior  of  a  symmetric  diblock 
copolymer  making  up  10  DPD  beads  and  consisting  of  5  A 
and  5  B  beads,  A5B5,  with  aAB  =  40  (which  corresponds  to 
Zab  —  2),  Ks  =  4,  and  r0  =  0.  Groot  and  Madden20  have  shown 
that  these  values  of  Ks  and  r0  provide  a  reasonable  description  of 
the  structural  properties  of  bulk  diblock  copolymer  melts  at 
p  =  3.  The  lengths  of  the  simulation  box  in  the  x  and  y  directions 
were  set  to  Lx  =  Lv  =  25,  and  the  slit  width  W  was  varied 
in  the  z  direction  from  2  to  20.  In  some  cases,  we  also  employed 
Lx  =  25  and  Ly  =  20  to  investigate  finite-size  and  commensurate 
effects  on  the  self-assembly,  typically  for  cases  near  or  at  the  phase 
boundaries  between  different  types  of  nanostructures.  Numerical 
integration  of  the  DPD  equations  of  motion  (eq  10)  was 
performed  using  the  modified  velocity  Verlet  algorithm.15  To 
save  CPU  time,  we  implemented  the  cell  list  method  for  the  energy 
and  force  calculations2122  Typically,  production  runs  spanned 
600  000  time  steps,  but  for  some  simulations  near  or  at  the  phase 
boundaries  between  different  types  of  nanostructures,  as  many  as 
2  000  000  time  steps  were  used. 

Simulations  of  the  spontaneous  formation  of  diblock 
copolymer  nanostructures  require  us  to  employ  a  convenient 
simulation  strategy  that  reasonably  guarantees  that  the  most 
stable  structure  will  be  found.  In  this  work,  we  use  the  following 
strategy  based  on  slow  cooling  of  the  diblock  copolymer  self- 
assembly.19  The  values  of  aAB  were  slowly  increased  from  aAB  = 
25  to  the  target  value  of  40,  whereby  because  ^ab  ^  1  /T23  this 
approach  corresponds  to  a  slow  decrease  in  T.  Each  increase  in 
aAB  was  followed  by  a  DPD  simulation  run  that  allowed  the 
system  to  relax  to  an  equilibrated  state.  The  increase  in  aAB  was 
typically  performed  with  steps  equal  to  unity  and  with  a  relaxa¬ 
tion  period  of  60  000  time  steps. 

In  addition  to  evaluating  thermodynamic  properties  such  as  the 
interparticle  energy,  the  pressure  tensor  (obtained  using  the 
Irvin— Kirkwood  definition),  the  temperature,  and  A-  and 
B-bead  density  profiles  across  the  slit,  we  also  evaluated  the 
following  structural  properties  of  the  diblock  copolymers  to 
deduce  their  conformational  behavior  in  confinement. 

The  mean  square  end-to-end  distance 
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and  the  mean  square  radius  of  gyration 
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are  typically  used  to  define  the  characteristic  size  of  a  polymer 
chain.  In  eqs  1 3  and  14,  TV  is  the  number  of  diblock  copolymers,  r,-  ; 
is  the  position  of  bead  /  belonging  to  diblock  copolymer  i,  rcom is 
the  position  of  the  center-of-mass  of  diblock  copolymer  i,  and 
(.)  denotes  an  ensemble  average.  Polymer  chains  confined  in 
planar  slits  can  be  further  characterized  by  the  mean  square 
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Clarendon  Press:  Oxford,  U.K.,  1987. 

(23)  Flory,  P.  J.  Principles  of  Polymer  Chemistry,  3rd  ed.;  Cornell  University  Press: 
New  York,  1953. 


radius-of-gyration  parallel  and  perpendicular  to  the  wall,  Rg j2 
and  Rg,_L,  respectively,  which  are  given  by  the  expressions 
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where  x,  y,  and  z  are  the  Cartesian  coordinates  of  r.  Note  that 
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Finally,  the  nematic  order  parameter 

3 

Rl,  d  —  2  (2|ti;ix) 


(16) 


is  a  convenient  measure  of  orientational  order,"4  which  we 
calculated  in  this  study.  In  eq  16,  Amax  is  the  largest  eigenvalue 
of  the  Saupe  tensor 
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where  e,-  is  the  unit  vector  along  the  end-to-end  direction  of 
diblock  copolymer  i  and  <5^  is  the  Kronecker  delta.  P2,a  is  zero  in 
a  completely  orientationally  disordered  state  and  is  unity  if  the 
diblock  copolymers  are  perfectly  aligned. 


3.  Results  and  Discussion 

We  carried  out  DPD  simulations  of  the  self-assembly  of  the 
A5B5  copolymers  in  planar  slits  with  and  without  nanopatterns. 
Knowledge  of  the  self-assembly  in  slits  without  nanopatterns 
provides  a  means  of  rationalizing  the  self-assembly  in  the  nano- 
patterned  slits.  We  were  particularly  interested  in  determining 
whether  the  self-assembly  behavior  in  the  slits  with  no  nanopat¬ 
terns  could  be  superimposed  to  predict  the  self-assembly  behavior 
accurately  in  the  various  nanopattern  scenarios.  Such  predictive 
capability  would  be  a  powerful  tool  in  the  rational  systematic 
design  of  nanostructures  and  nanostructure  templates. 

In  this  study,  we  varied  the  slit  width  W from  2  to  20  with  a  unit 
step  and  varied  the  DPD  repulsion  parameters  aAW]  and  aBH,2 
typically  from  the  value  of  aAA  =  aBB  to  the  value  of  aAB  (i.e.,  from 
25  to  40).  Ratios  of  AaAn.JAaAB  and  AaBM,2/AflAB,  where  AaAv/i  = 
Qak.,  “  aAA,  A aBM,2  =  oBh,2  -  aAA,  and  AaAB  =  uAb  “  «aa,  were 
used  to  characterize  the  phobicity  of  the  A  and  B  beads  toward  the 
wr  and  W2-wall  beads,  respectively.  The  phobicity  increases  with 
increasing  values  of  AaAwJAaAB  or  analogously  AflBH,2/AaAB. 
Note  that  AaAu,2  =  Aas,,^  =  0  for  our  chosen  parameter  set 
given  in  section  2.1. 

3.1.  Planar  Slits  without  Nanopatterns.  We  first  studied 
the  self-assembly  of  A5B5  copolymers  confined  between  planar 
walls  where  the  opposing  walls  are  (a)  both  composed  of  uq-wall 
beads  (referred  to  hereafter  as  W\  walls)  and  (b)  composed  of 
different  types  of  wall  beads,  w’i  wall  and  w’2  wall,  respectively 
(referred  to  hereafter  as  wq  and  w2  walls). 

3.1.1.  Planar  Slits  with  w /  Walls.  We  studied  systems  with 
repulsion  parameters  AaAn.JAaAB  =  {0,  73,  2/s,  1, 4/3}  where  for 
all  of  these  systems,  except  the  first  listed,  B  beads  have  a 
preferential  affinity  for  uq-wall  beads  because  Aas^  =  0.  Figure  2 
displays  an  aAw  —  W  phase  diagram  for  the  A5B5  system,  where 
the  particular  phases  were  determined  by  an  analysis  of  both  the 
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Figure  2.  Phase  diagram  of  the  A5B5  system  confined  in  planar  slits  with  w  i  walls,  together  with  examples  of  simulation  configurations  for  the 
particular  lamellar  phases.  Wis  the  slit  width,  and  AaAn.JAaAB  is  a  measure  of  the  phobicity  of  the  A  beads  toward  the  w^-wall  beads;  B  beads 
have  preferentially  affinity  toward  the  H’i-wall  beads  (Aubw/Auab  =  0).  Symbols  represent  simulation  results,  and  solid  lines  indicate 
approximate  microphase  separation  boundaries.  The  symbols  refer  to  a  perpendicular  lamellar  phase  (LjJ  and  parallel  lamellar  phases  with 
v  A  lamellae  (Ln  „).  Green  and  gold  spheres  represent  A  and  B  beads,  respectively,  and  gray  spheres  represent  wall  beads. 


isosurfaces  and  the  A-  and  B-bead  density  profiles  across  the 
slit,  together  with  a  visual  inspection  of  the  equilibrium 
configurations.  For  narrow  slit  widths  up  to  4  or  5,  the  system 
forms  a  perpendicular  lamellar  phase  (denoted  as  L  x)  regardless  of 
the  interaction  strength  between  the  chains  and  the  walls  for  steric 
reasons  (i.e.,  the  slit  widths  are  simply  too  narrow  to  accommodate 
parallel  lamellae).  Figure  2  also  shows  that  when  no  preferential 
affinity  of  the  A  or  B  beads  toward  the  wall  beads  exists 
(i.e.,  AaAwJAaAB  =  0)  the  system  self-assembles  into  a  perpendi¬ 
cular  lamellar  phase  for  all  slit  widths.  For  these  cases,  perpendi¬ 
cular  alignment  is  preferred  over  parallel  alignment  because  the 
wall-diblock  copolymer  adsorption  energy  does  not  play  a  sig¬ 
nificant  role;  moreover,  the  perpendicular  alignment  allows  the 
system  to  adopt  its  preferred  bulk  lamellae  spacing.  Figure  2 
further  shows  that  with  increasing  AaAwJAaAB  (i.e.,  phobicity  of 
the  A  beads  toward  the  tv  i  -wall  beads)  the  systems  at  W  >  5  prefer 
to  form  parallel  lamellar  phases  spontaneously  with  v  lamellae 
(denoted  as  L||v)  as  a  result  of  the  more  favorable  wall— diblock 
copolymer  adsorption  energy  for  parallel  alignment  with  respect  to 


perpendicular  alignment,  which  is  a  consequence  of  the  preferential 
adsorption  of  the  B  beads.  (Note  that  throughout  this  article 
v  refers  to  the  number  of  A  lamellae  unless  indicated  otherwise;  for 
example,  in  Figure  la,  v  =  2.)  The  lamellae  adopts  a  perpendicular 
alignment  only  when  the  system  becomes  very  structurally 
frustrated  between  the  L||iV  and  LnjV+1  orientations.  In  such  cases, 
the  slit  width  does  not  correspond  to  either  the  v  or  the  v  +  1 
lamellae  spacing.  For  strong  phobicity  between  the  A  beads  and 
wall  (AaAn,JAaAB  >  0.85),  Lx is  formed  only  between  Fi.i  and  La  2. 

Figures  3  and  4  show  R2,  Rg 2,  and  P2, a  as  a  function  of  W  for 
two  particular  values  of  phobicity,  AaAn,JAaAB  =  0  and  AaAwJ 
AaAB  =  1,  which  provide  further  insight  into  the  behavior  of  the 

diblock  copolymer  self-assembly.  For  AflAu.yAflAB  =  0,  where  the 

system  forms  only  perpendicular  lamellae,  the  diblock  copoly¬ 
mers  become  increasingly  elongated  in  the  perpendicular  direc¬ 
tion  and  the  system  loses  orientational  order  (i.e.,  P2,a  decreases) 
as  W  increases  up  to  about  6.  Because  the  chains  in  this  case  have 
no  preferential  interactions  with  the  walls,  this  behavior  can  be 
attributed  only  to  a  release  of  the  steric  constraint  imposed  by  the 
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Figure  3.  (a)  Mean  square  end-to-end  distance,  /?e2,  (b)  mean  square  radius-of-gyration,  f?g2,  and  (c)  nematic  order  parameter,  P2, d,  as  a 
function  of  the  slit  width,  W,  for  the  A5B5  copolymers  confined  in  planar  slits  with  W|  walls,  where  Aa^JAa^B  =  Aub^/A^ab  =  0.  Symbols 
represent  simulation  results,  where  dotted  lines  are  drawn  as  a  guide  to  the  eye. 


confining  walls.  Above  W~  6,  the  size  of  the  copolymers  and  the 
system  orientational  order  is  only  marginally  affected  by  W. 

Very  different  conformational  behavior  for  the  diblock  copo¬ 
lymers  is  seen  for  AaAlVi/AaAB  =  1  in  Figure  4,  where  the  system 
self-assembles  into  both  parallel  and  perpendicular  alignments 
depending  on  W.  For  W  <  6,  where  the  systems  exhibit 
perpendicular  alignment,  R2,  Pg2,  and  P2,d  exhibit  similar  trends 
as  in  the  case  of  AaAwJAaAB  =  0  (i.e.,  the  diblock  copolymers 
elongate  in  the  perpendicular  direction  and  the  system  orienta¬ 
tional  order  decreases).  For  W  >  6,  the  transformation  into  a 
particular  parallel  phase  is  accompanied  by  a  dramatic  elongation 
of  the  diblock  copolymers  in  the  perpendicular  direction,  as 
evidenced  by  the  increase  in  R2,  Rg  j_2,  and  P2,d  with  increasing 


W  within  each  Lii,,-  phase.  As  the  slit  width  is  increased,  structural 
frustration  of  the  diblock  copolymers  in  the  parallel  lamellae  is 
relieved  by  either  the  formation  of  an  additional  parallel  lamellae 
or,  in  the  case  where  W  =*  10,  a  reorientation  of  the  parallel 
lamellae  to  perpendicular  lamellae.  For  either  of  these  occur¬ 
rences,  the  diblock  copolymers  have  collapsed  upon  themselves  to 
minimize  the  free  energy  while  subsequently  losing  orientational 
order.  Note  that  for  W  >  6,  Pgj[ 2  adopts  distinct  values  for  Lj_  and 
Lj|  r  and  they  do  not  depend  on  W  within  the  perpendicular 
and  parallel  lamellar  regions.  For  the  other  nonzero  values  of 
AaA„,  /AaAB  shown  in  Figure  2,  we  observed  similar  dramatic 
differences  in  the  conformational  behavior  when  compared  to 
AaAvV|/A<7AB  =  0,  but  they  are  not  shown  here. 
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Figure  4.  (a)  Mean  square  end-to-end  distance,  R,.2.  (b)  mean  square  radius  of  gyration,  Rg2,  and  (c)  nematic  order  parameter,  P2,d,  as  a 
function  of  the  slit  width,  W ,  for  the  A5B5  copolymers  confined  in  planar  slits  with  uq  walls,  where  AaAwJAaAB  =  1  and  AaB„,JAaAB  =  0. 
Symbols  represent  simulation  results  where  dotted  lines  are  drawn  as  a  guide  to  the  eye  and  vertical  dashed  lines  mark  approximate 
microphase  separation  boundaries.  Symbols  for  particular  phases  are  the  same  as  in  Figure  2. 


The  conformational  behavior  exhibited  in  Figures  3  and  4 
correlates  with  the  dependence  of  the  lamellae  width,  WLAM,  on 
W  as  shown  in  Figure  5,  where  the  values  of  WLAM  were 
determined  by  Voronoi  tessellation.25  Note  that  for  W  <  6, 
perpendicular  lamellae  are  quite  irregular  and  therefore  the  value 
of  WLAM  cannot  be  reliably  determined.  As  expected  for  A aAwJ 
Aoab  =  0,  the  lamellar  width  does  not  depend  on  W  and  the 
values  of  IFLam  correspond  to  the  natural  lamellar  width  in  the 


(25)  Arya,  S.;  Malamatos,  T.;  Mount,  D.  M.  Space- Efficient  Approximate 
Voronoi  Diagrams, ;  Proceedings  of  the  34th  ACM  Symposium  on  Theory  of 
Computing;  STOC,  Monreal,  Canada,  2002;  pp  721—730. 


bulk.  For  AaAwJAaAB  =  1,  the  lamellar  widths  increase  with 
increasing  W  within  each  parallel  lamellar  region  of  the  phase 
diagram  coinciding  with  the  elongation  of  the  chains  in  the 
perpendicular  direction  and  the  increase  in  orientational  order 
exhibited  in  Figures  3  and  4.  For  the  parallel  lamellae  phases  only, 
Figure  5b  also  displays  the  thickness  of  the  adsorbed  B-bead 
layer,  WL,  where  the  values  of  1FL  were  evaluated  to  be  W- L  = 
j2[W  —  (2v  —  \)W\  AM].  Similar  to  the  behavior  found  for  the 
lamella  thickness,  we  see  that  with  increasing  slit  width  the 
adsorbed  layer  thickness  increases  within  each  parallel  lamellar 
region,  however  in  a  less  dramatic  manner.  In  addition,  the 
adsorbed  layer  thickness  becomes  less  dependent  on  the  slit  width 


Langmuir  XXXX,  XXX (XX),  XXX-XXX 


DOI:  10.1 02 1  /la903200j  G 


Downloaded  by  Martin  Lisal  on  October  19,  2009  I  http://pubs.acs.org 
Publication  Date  (Web):  October  19,  2009  I  doi:  10.1021/la903200j 


Article 


Petrus  et  al. 


1-11,1  L|  ^11,2  ^11,3 


W 

Figure  5.  Lamella  width,  WLAM,  and  thickness  of  the  adsorbed  layer,  WL,  as  a  function  of  the  slit  width,  W ,  for  the  A5B5  copolymers  confined 
in  planar  slits  with  w i  walls,  where  (a)  AciAwJAaAB  =  AaBwJAaAB  =  0  and  (b)  AaAwJAaAB  =  1  and  AaBwJAaAB  =  0.  Symbols  represent 
simulation  results,  where  dotted  lines  are  drawn  as  a  guide  to  the  eye  and  vertical  dashed  lines  mark  approximate  microphase  separation 
boundaries.  Symbols  for  particular  phases  are  the  same  as  in  Figure  2. 


when  more  lamellae  are  present  in  the  system.  This  behavior  can 
be  attributed  to  the  decreasing  influence  of  the  wall  interactions 
because  as  the  slit  width  increases,  the  system  is  approaching  bulk- 
phase  behavior. 

3.1.2.  Planar  Slits  with  w ,  and  w2  Walls.  Figure  6  presents  an 
flBivj-  W  phase  diagram  for  the  confined  A5B5  copolymers,  where 
the  A  beads  interact  with  the  w\  walls  by  AaAlV|/AaAB  =  1  and  the 
B —w2  repulsion  interaction  is  varied  as  AaBwJAaAB  =  {0,  '/2,  1}. 
Here,  A  and  B  beads  have  preferential  affinity  for  the  w2-  and 
M’]-wall  beads,  respectively,  because  AaAll,  =  A aBWi  =  0.  As  in  the 
previous  case  where  the  opposing  walls  were  identical,  for  very 
narrow  slit  widths  of  up  to  2  or  3  the  system  self-assembles  into 
perpendicular  alignment  regardless  of  the  strength  of  the  interac¬ 
tion  between  the  diblock  copolymers  and  the  walls.  For  narrow  slit 
widths  between  3(4)  and  5,  the  system  forms  only  A  and  B  adsorbed 
layers,  which  we  denote  as  an  adsorbed-layer  phase,  L.  As 
another  similarity  to  the  previous  case,  a  variety  of  parallel 
alignments  are  formed  at  larger  slit  widths.  However,  the  system 


displays  two  distinguishable  phases  with  parallel  alignment:  parallel 
lamellar  phases  with  B-adsorbed  layers  on  both  confining  walls 
(as  in  the  previous  case)  and  parallel  lamellar  phases  with 
A-  and  B-layers  adsorbed  onto  the  vv2  and  W\  walls,  respectively, 
with  the  latter  phases  termed  adsorbed  layer— parallel  lamellar 
phases  and  denoted  as  L— b\.v.  In  our  simulations,  L— Lii,„  phases 
exist  for  AaBwJAaAB  <  1,  where  B  beads  can  also  adsorb  onto  the 
w2  wall  because  of  the  lower  B-wall  phobicity  with  respect  to  the 
A-wall  phobicity  (i.e.,  due  to  AaBvl,2  <  AaAWi). For  AaBwJAaAB  =  0, 
L  and  L— L||  „  phases  always  precede  Hi.,,  and  Lh,v+1  phases, 
respectively.  However,  for  AaBwJAaAB  =  0.5,  we  see  the  disap¬ 
pearance  of  L|i,2  and  L|i,3  phases,  along  with  the  formation  of  Lx 
between  L  and  Luj,  and  a  mixed  lamellar  phase  (Lx+u)  between 
L —  L||  i  and  L — Lu  <2.  For  all  slit  widths,  the  system  with  AaBMS/AaAB  = 
AaAwJAaAB  =  1  does  not  exhibit  Lii.v  (as  expected  from  the 
symmetry  in  the  diblock  copolymer— wall  interactions)  and  displays 
a  succession  of  Lx/L/Lx/L— Ln^/L—  Lu  j2/L —  Lh3  phases.  Note  that 
we  have  verified  the  stability  of  the  perpendicularly  aligned  phases 
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Figure  6.  Phase  diagram  of  the  A5B5  system  confined  in  planar  slits  with  m’!  and  vy2  walls,  where  AaAwJAaAB  =  1  and  AaBwJ  AaAB  varies. 
Also  shown  are  examples  of  simulation  configurations  for  the  layer,  layer— parallel  lamellar,  and  mixed  lamellar  phases;  simulation 
configurations  for  the  perpendicular  and  parallel  lamellar  phases  are  analogous  to  those  shown  in  Figure  2.  W  is  the  slit  width,  and 
AaAwJAaAB  and  A<7Bm,2/Aoab  are  measures  of  the  phobicity  of  the  A  and  B  beads  toward  the  u’r  and  it’2-wall  beads,  respectively.  Symbols 
represent  simulation  results,  and  solid  lines  mark  approximate  phase  boundaries.  Symbols  refer  to  a  perpendicular  lamellar  phase  (LJ,  a 
mixed  lamellar  phase  (Lx -Ml),  parallel  lamellar  phases  with  v  A  lamellae  (Ll.v),  a  layer  phase  (L),  and  layer— parallel  lamellar  phases  with  v  A 
lamellae  (L— L|i  v).  Green  and  gold  spheres  represent  A  and  B  beads,  respectively,  and  gray  and  pink  spheres  represent  w i-  and  vt>2-wall  beads, 
respectively. 


for  W  >  5  and  the  mixed  lamellar  phase  Lx+ii  by  performing  DPD 
simulations  with  box  lengths  of  Lx  =  25  and  Lv  =  20. 

As  an  example  of  the  conformational  behavior  of  the  systems 
shown  in  Figure  6,  Rc2,  Rg2,  and  P2a  are  shown  in  Figure  7  for 
AaBwJAaAS  =  0.  Similar  to  the  behavior  exhibited  in  planar  slits 
with  W|  walls  when  AaAwJAaAB  =  1  and  AaBlVi/AaAB  =  0,  the 
succession  of  L/L|i,i/L— L||,i/Lii,2/L—  L|i,2/L|i,3  phases  is  accompa¬ 
nied  by  the  elongation  of  the  diblock  copolymers  in  the  perpen¬ 
dicular  direction  and  an  increase  in  the  system  orientational  order 
within  each  phase  for  W  >  3.  The  stretching  and  increasing 
orientational  order  is  then  relieved  by  the  transformation  of  one 
phase  into  another  as  W  increases.  Similar  trends  in  the  con¬ 
formational  behavior  as  depicted  by  R2,  Rg 2,  and  P2,d  are  also 
exhibited  for  AflBu,2/AflAB  >  0  (not  shown  here). 

3.2.  Nanopatterned  Planar  Slits.  Next,  we  study  the  self- 
assembly  of  A5B5  copolymers  confined  in  nanopatterned  planar 
slits  for  three  mutual  locations  of  the  nanopattems  and  for  two 


differently  sized  nanopatterned  strips:  LNP  =  12.5  and  6.25,  as 
depicted  in  Figure  8.  In  Figure  8,  opposing  walls  have  strips 
composed  of  wr  and  w2- wall  beads  that  are  vertically  (a)  aligned, 
(b)  staggered,  and  (c)  partially  staggered.  After  setting  the 
A-bead  repulsion  parameters  to  AuAM,i/AaAB  =  1  and  A aAw2  = 
AaHlV|  =  0,  we  studied  systems  for  B-bead  repulsion  parameters 
AflBw,/AflAB  =  {0,  V2,  !}■ 

In  Figure  8,  the  nanopatterns  on  opposing  walls  can  be 
partitioned  such  that  they  correspond  to  the  planar  surfaces 
without  nanopattems  investigated  previously  and  are  indicated 
as  such  by  the  red  dashed  lines.  Hence,  it  is  reasonable  to  attempt 
to  rationalize  the  self-assembly  of  the  A5B5  copolymers  in  the 
nanopatterned  planar  slits  on  the  basis  of  the  behavior  observed 
in  the  planar  slits  without  nanopatterns.  In  particular,  we 
investigate  the  applicability  of  superimposing  the  observed  nano¬ 
structures  in  planar  slits  without  nanopatterns  to  predict 
the  nanostructures  formed  when  nanopatterns  are  present. 
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Figure  7.  (a)  Mean  square  end-to-end  distance,  Re2,  (b)  mean  square  radius  of  gyration,  7?g2,  and  (c)  nematic  order  parameter,  P2,d>  as  a 
function  of  the  slit  width,  W ,  for  the  A5B5  copolymers  confined  in  planar  slits  with  w  i  and  w2  walls  and  with  AuawJAqab  =  1  and  AubwJ  Agab 
=  0.  Symbols  represent  simulation  results  where  dotted  lines  are  drawn  as  a  guide  to  the  eye  and  vertical  dashed  lines  mark  approximate 
microphase  separation  boundaries.  Symbols  for  particular  phases  are  the  same  as  in  Figure  6. 


Corresponding  to  the  partitions  indicated  in  Figure  8,  the  super¬ 
position  description  is  satisfied  if  the  nanostructure  in  the  nano- 
patterned  slit  can  be  approximated  by  a  collective  sequence  of 
previously  determined  nanostructures.  Our  attempt  to  depict  a 
phase  diagram  via  such  a  superposition  description  is  performed  in 
the  following  manner.  For  each  slit  width,  B-w2  phobicity  pair  in 
the  phase  diagram,  we  first  provide  a  schematic  indicating  the  type 
of  nanostructure  observed  in  the  planar  slits  without  nanopatterns 
that  corresponds  to  each  of  the  partitions.  If  a  nanostructure  for  a 
particular  W—  AaBw2  pair  corresponds  to  simply  superimposing 
the  nanostructures  observed  in  planar  slits  without  nanopatterns, 


then  we  provide  an  equilibrium  configuration  snapshot  (exceptions 
are  for  previously  shown  configurations  such  as  Lj^  or  fi„). 
Otherwise,  a  cross  is  placed  over  the  schematic  of  the 
nanostructure(s)  corresponding  to  the  case  without  nanopatterns, 
whereas  next  to  this  schematic  we  indicate  the  actual  nanostruc¬ 
ture  that  is  formed.  In  other  words,  in  the  phase  diagram,  the 
schematics  of  nanostructures  that  have  a  cross  through  them  do 
not  follow  from  superimposing  the  analogous  nanostructures 
formed  by  the  self-assembly  in  planar  slits  without  nanopatterns. 

As  a  further  explanation  of  constructing  the  schematics  that  are 
used  to  depict  the  phase  behavior,  Figure  9  provides  three 


J  DOI:  10.1021/la903200j 


Langmuir  XXXX,  XXX (XX),  XXX  XXX 


Downloaded  by  Martin  Lisal  on  October  19,  2009  I  http://pubs.acs.org 
Publication  Date  (Web):  October  19,  2009  I  doi:  10.1021/la903200j 


Petrus  et  a!. 


Article 


Figure  8.  Nanopatterned  planar  slits,  where  the  nanopatterns  on 
opposing  surfaces  are  vertically  (a)  aligned,  (b)  staggered,  and  (c) 
partially  staggered.  Left  and  right  columns  correspond  to  the 
nanopatterned  planar  slits  with  wide  (Z.NP  =  12.5)  and  narrow 
(Z,NP  =  6.25)  strips,  respectively.  Red  dashed  lines  indicate  virtual 
partitioning  of  the  nanopatterned  slits  into  portions  corresponding 
to  the  planar  slits  without  nanopatterns  (i.e.,  planar  slits  with  vi’i 
walls  and  planar  slits  with  w\  and  w2  walls).  Gray  and  pink  spheres 
represent  u'r  and  wj-wall  beads,  respectively. 

examples  of  some  nanostructures  that  are  observed  in  planar  slits 
with  wide,  vertically  aligned  nanopatternes.  In  the  first  example 
where  W  =  7  (Figure  9a),  the  schematic  indicates  the  formation 
of  L,Ui  and  L±  in  the  w \  and  w2  portions  of  the  nanopatterned 
planar  slit,  respectively.  The  total  phase  behavior  in  the  nano- 
patterned  planar  slit  can  be  predicted  by  superimposing  the  Li.i 
and  L±  nanostructures  observed  in  planar  slits  without  nanopat¬ 
terns,  more  precisely,  in  planar  slits  of  W  =  7  that  exhibit  L±  and 
Li,i  phases  for  0  and  1  block-wall  phobicities,  respectively, 
as  shown  in  Figure  2.  In  the  second  example  where  W  = 
20  (Figure  9b),  the  schematic  indicates  that  the  superposition 
description  is  not  valid  because  instead  of  the  formation  of  L|i,3 
and  L±  phases  in  the  wt  and  w2  portions  of  the  nanopatterned 
planar  slit,  respectively,  the  system  exhibits  Li,3  behavior  only. 
Therefore,  in  this  particular  case,  the  nanopatterned  system  only 
partially  follows  the  superposition  description,  where  a  cross 
overlays  the  part  of  the  schematic  corresponding  to  L±.  In  the 
third  example  where  IV  =  1 1  (Figure  9c),  the  behavior  observed 
in  the  planar  slits  without  nanopatterns  would  suggest  self- 
assembly  into  an  L|i,2  phase  and  an  inverse  L|i,2  phase  (indicated 
by  Lii,2_1)  in  the  ID  and  w2  portions  of  the  nanopatterned  planar 
slits,  respectively,  but  instead  the  system  self-assembles  into  only 
an  L±  phase.  (In  an  Li,v  1  phase,  the  roles  of  the  A  and  B  beads 
are  interchanged,  where  now  the  A  beads  exhibit  preferential 
affinity  toward  the  walls  and  lamellae  are  formed  by  B  beads  (and 
correspondingly  v  refers  to  the  number  of  B  lamellae).)  Here,  the 
system  does  not  follow  the  superposition  description  in  either  the 
vv  |  or  the  w2  portion  of  the  nanopatterned  planar  slit,  thus  both 
parts  of  the  schematic  are  crossed  out. 

As  performed  for  planar  slits  without  nanopatterns,  in  addition 
to  DPD  simulations  with  box  lengths  of  Lx  =  Ly  =  25,  we  also 
investigated  finite-size  and  commensurate  effects  by  performing 
additional  DPD  simulations  with  Lx  =  25  and  Ly  =  20.  In  most 
cases,  these  simulations  confirmed  our  findings  from  the  original- 
size  simulations,  except  for  a  few  cases  near  the  microphase- 
separation  boundaries.  In  such  cases,  we  indicate  both  self- 
assembled  nanostructures  in  the  phase  diagrams.  A  determination 
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Figure  9.  Examples  of  the  approach  used  to  construct  the  sche¬ 
matic  phase  diagrams  of  A5B5  copolymers  confined  in  planar  slits 
with  wide,  vertically  aligned  nanopatterns.  In  portions  a  and  b,  B 
beads  have  preferential  affinity  toward  vv  i  strips  and  neither  A  nor 
B  beads  have  preferential  affinity  toward  w2  strips.  In  portion  c,  A 
and  B  beads  have  preferential  affinity  toward  w2  and  w i  strips, 
respectively.  The  left  column  shows  examples  of  the  schematics 
that  are  reported  in  the  phase  diagram,  whereas  the  right  column 
shows  the  observed  phase  behavior  via  a  simulation  configuration. 
Each  schematic  consists  of  two  parts  corresponding  to  w\  and  w2 
portions  of  the  nanopatterned  slits,  where  for  each  partition  the 
lamellar  phase  found  in  the  corresponding  planar  slit  without 
nanopatterns  is  displayed.  Crosses  through  a  partition  of  the 
schematic  indicate  that  this  phase  is  not  present  in  the  actual 
nanopatterned  planar  slit  and  thus  a  superposition  description  is 
not  valid.  Green  and  gold  spheres  represent  A  and  B  beads, 
respectively,  and  gray  and  pink  spheres  represent  w\  and  w2  wall 
beads,  respectively. 

of  the  true  nanostructures  very  near  these  microphase  separation 
boundaries  would  require  a  systematic  analysis  of  finite-size  and 
commensurate  effects  accompanied  by  free-energy  calculations,26 
which  is  beyond  the  scope  of  this  work. 

3.2.1.  Planar  Slits  with  Wide,  Vertically-Aligned  Nanopat¬ 
terns.  Figure  10  presents  the  phase  behavior  of  A5B5  copolymers 
confined  in  planar  slits  with  wide  (Z,NP  =  12.5),  vertically  aligned 
nanopatterns  in  a  schematic  manner  for  the  various  pairs  of  W 
and  AflB„,2/AflAB  studied.  Figure  1 1  provides  simulation  config- 


(26)  Martinez-Veracoeche,  F.  J.;  Escobedo,  F.  A.  J.  Chem.  Phys.  2006,  125, 
104907. 
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Figure  10.  Self-assembly  of  the  A5B5  systems  confined  in  nano- 
patterned  planar  slits,  where  strips  on  opposing  surfaces  are 
vertically  aligned  with  size  Z,NP  =  12.5.  W  is  the  slit  width,  and 
AaB>v2/AflAB  is  a  measure  of  the  phobicity  of  the  B  blocks  toward 
the  h’2  walls.  Schematics  indicate  the  lamellar  alignment  in  the 
planar  slits  without  nanopattems  (i.e.,  in  the  planar  slits  with  vt’i 
walls),  w i  and  u'2  walls  are  represented  by  black  and  blue  lines, 
respectively.  Crosses  indicate  cases  where  nanostructures  did  not 
result  from  the  superposition  of  the  lamellar  alignment  in  the 
underlying  planar  slits  without  nanopatterns.  Letters  a— g  repre¬ 
sent  simulation  configurations  for  some  select  nanostructures 
formed  in  the  nanopatterned  planar  slits,  which  are  shown  in 
Figure  1 1 .  Other  legend  definitions  are  Lx,  a  perpendicular  lamellar 
phase,  and  Ln  3,  a  parallel  lamellar  phase  with  three  A  lamellae.  For 
cases  in  which  different  nanostructures  were  observed  for  differ¬ 
ently  sized  simulation  cells,  both  nanostructures  are  shown,  sepa¬ 
rated  by  a  comma. 

urations  for  some  of  the  nanostructures  formed,  which  corres¬ 
pond  to  the  letters  given  in  Figure  10.  A  few  interesting  features  of 
Figure  10  are  worth  noting.  Overall,  the  resulting  nanostructures 
follow  from  the  superposition  of  the  lamellar  alignment  observed 
in  the  planar  slits  without  nanopatterns  (i.e.,  in  the  planar  slits 
with  w’i  walls),  with  an  exception  for  wide  slit  widths,  W  >  15.  For 
narrow  slit  widths  of  up  to  5,  the  system  follows  superposition 
behavior,  where  Lx  phases  are  formed  regardless  of  the  value  of 
AflBM,2/Af/AB.  Considering  all  slit  widths  for  the  A«BH,2/AaAB  =  0 
case,  we  observed  the  formation  of  nanostructures  given  by  the 
superposition  behavior  of  the  L|i,v  and  Lx  phases  in  all  cases 
(denoted  as  U\,v  +  Lj_,  where  configuration  snapshots  are  shown 
in  Figure  11a— c),  with  the  exception  of  the  relaxation  of 
the  system  into  Lx  phases  between  L|[,i  +  Lx  and  Li,2  +  Lx  and 
between  Ln  2  +  Lx  and  Ln,3  +  Lj_  phases.  The  observed  interface 
morphology  between  parallel  and  perpendicular  lamellar 
domains  is  known  as  Schrek’s  surface  twist  boundary.  Schrek’s 


Figure  11.  Examples  of  simulation  configurations  corresponding 
to  Figure  10.  (a)  L,i  +  Lx,  (b)  U\,i  +  Lx,  (c)  Lii,3  +  Lx,  (d)  Ln,,  + 
Ll  i  ,  (e)  Cl  ,2  +  Cl  ,2  1 ,  (f)  tilt-like  lamellar  phase,  and  (g)  chevron¬ 
like  lamellar  phase.  Shown  in  the  inset  of  panel  a  are  isosurfaces  for 
the  A  beads  together  with  Schrek’s  surface  twist  boundary  indicted 
by  a  red  oval.  Green  and  gold  spheres  represent  A  and  B  beads, 
respectively,  and  gray  and  pink  spheres  represent  w\  and  w2  wall 
beads,  respectively.  B  beads  are  not  shown  in  panels  d  and  e  for 
visual  clarity. 

surface  is  generated  by  smoothly  joining  two  sets  of  planes  that 
are  aligned  normal  with  respect  to  each  other.27'28  (Schrek’s 
surface  twist  boundary  for  the  L|i,i  +  Lx  phase  is  seen  in  the  inset 
of  Figure  11a.)  The  interface  morphology  of  Schrek’s  surface 
minimizes  the  energy  penalty  at  the  boundary  interface  and 
thereby  minimizes  the  overall  energy  of  the  A5B5  copolymers  in 
the  L|i,v  +  L_l  phases.12  The  L|i,v  +  Lx  phases  prefer  Schrek’s 
surface  twist  boundary  over  a  T-junction  tilt  boundary,29  again 
because  of  the  lower  energy  penalty  at  the  boundary  interface.  For 
W  =  (19)20  when  A«Bm,2/Aaab  =  0,  the  influence  of  the  w2  walls  is 
now  so  weak  with  respect  to  the  other  interactions  that  the  system 
prefers  to  self-assemble  into  a  parallel  alignment,  closely  resem¬ 
bling  its  phase  behavior  in  the  planar  slits  without  nanopatterns. 

Overall,  for  the  cases  of  AaBwJAaAB  =  0.5  and  1,  once  again, 
superposition  behavior  is  followed  where  nanostructures  with 
alternating  parallel  and  inverse  parallel  lamellae  are  found 
(denoted  as  L|i,i  +  Ln  j-1  and  U\,2  +  Li,2  >  where  configuration 
snapshots  are  shown  in  Figure  1  ld,e,  respectively).  As  observed 
previously  in  the  AflBH,2/AaAB  =  0.0  case,  relaxation  of  the  system 
into  perpendicular  alignment  between  Li,i  +  Ln^-1  and  Li,2  + 
L|i,2_1  is  a  result  of  the  slit  width  not  being  commensurate  with  the 
bulk  lamellae  spacing.  As  also  observed  previously,  for  wide  slit 


(27)  Thomas,  E.  L.;  Anderson,  D.  M.;  Henkee,  C.  S.;  Hoffman,  D.  Nature  1988, 
334,  598. 

(28)  Gido,  S.  P.;  Gunther,  J.;  Thomas,  E.  L.;  Hoffman,  D.  Macromolecules  1993, 
26,  4506. 

(29)  Gido,  S.  P.;  Thomas,  E.  L.  Macromolecules  1994,  27,  6137. 
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widths  the  superposition  behavior  fails.  For  W  =  18,  19,  and 
20  when  AobhJAuab  =  0.5,  the  systems  self-assemble  into  tilt-like 
lamellar  phases  as  shown  in  Figure  1  If.  For  IF  between  16  to  20 
when  AflBM,2/AaAB  =  1,  we  observed  chevron-like  lamellar  phases 
(Figure  1  lg)  instead  of  the  expected  Lij  +  Li,3  •  The  simulation 
results  in  these  wider  slits  indicate  a  weakened  influence  of  the 
nanopatterns  on  the  self-assembly  because  the  tilt-like  and 
chevron-like  lamellae  are  similar  to  the  Li,v  phases  formed  in 
planar  slits  without  nanopatterns. 

Additional  DPD  simulations  with  box  lengths  of  Lx  =  25  and 
Ly  =  20  confirmed  the  self-assembled  nanostructures  obtained  by 
the  simulations  with  Lx  =  Ly  =  25  except  in  the  cases  of  [  W , 
AuBh,7Auab]  —  {[16,0],  [19,0],  [9,0.5]},  where  the  simulations 
relaxed  to  Li,2  +  Lx,  Li  ,3  +  Lx,  and  Lx,  respectively.  Nano¬ 
structures  determined  from  both  simulation  sizes  are  indicated  in 
Figure  10,  where  the  original  size  simulation  result  is  shown  first, 
separated  by  comma.  All  of  these  cases  occur  at  microphase 
separation  boundaries  between  two  different  types  of  nanostruc¬ 
tures,  which  explains  the  ambiguity  due  to  highly  pronounced 
finite-size  and  commensurate  effects. 

The  schematic  phase  diagram  and  some  select  nanostructures 
for  narrow  (Lnp  =  6.25),  vertically  aligned  nanopatterns  together 
with  a  detailed  discussion  are  left  for  the  Supporting  Information. 
Furthermore,  examples  of  conformational  behavior  in  terms  of 
the  nematic  order  parameter  for  the  A5B5  copolymers  confined  in 
the  vertically  aligned  nanopattern  slits  are  also  provided  in  the 
Supporting  Information. 

3.2.2.  Planar  Slits  with  Wide,  Vertically-Staggered  Nano¬ 
patterns.  The  schematic  phase  diagram  and  some  select 
nanostructures  for  wide  (Lnp  =  12.5),  vertically  staggered  nano¬ 
patterns  are  presented  in  Figures  12  and  13,  respectively.  Analo¬ 
gous  to  planar  slits  with  vertically  aligned,  wide  nanopatterns, 
generally  the  system  self-assembles  into  nanostructures  given  by 
the  superposition  behavior  of  lamellar  alignment  in  planar  slits 
without  nanopatterns  (i.e.,  in  the  planar  slits  with  vv |  and  w2 
walls),  with  some  exceptions  for  AflBu.2/AflAB  —  0.5.  We  observed 
the  formation  of  nanostructures  given  by  the  superposition  of  L 
and  inverted  L  (denoted  as  L  +  iL,  where  a  configuration 
snapshot  is  given  in  Figure  13a)  and  the  superposition  of  L— Li,„ 
and  inverted  L—  Lii.v  with  v  =  1  and  2  (denoted  as  L— Lii.v  + 
iL— Li,„,  where  configuration  snapshots  at  two  different  slit 
widths  are  given  in  Figure  13b,c).  (Inverted  adsorbed  layer  or 
inverted  adsorbed  layer— parallel  lamellar  phases  differ  from 
adsorbed  layer  or  adsorbed  layer— parallel  lamellar  phases  simply 
by  the  adsorbed  layer  positioning  with  respect  to  the  confining  w i 
and  u'2  walls.  Thus,  in  our  simulations,  the  adsorbed  layer  and 
adsorbed  layer— parallel  lamellar  phases  are  confined  by  the 
upper  w i  wall  and  lower  w2  wall,  and  the  positions  of  the  wt 
and  w2  walls  are  inverted  for  the  so-called  inverted  phases.) 
A  unique  nanostructure  was  observed  for  [W,  AaHwJAa^ B]  = 
[8,  0.5],  where  the  system  self-assembles  into  U-shape  lamellae 
shown  in  Figure  13d.  (This  U-shape  lamellar  phase  was  also 
confirmed  by  a  DPD  simulation  with  box  lengths  of  Lx  =  25  and 
Ly  =  20.)  The  U-shaped  lamellar  phase  is  evidently  a  transition 
structure  between  Lx  and  L— Li,„  +  iL— Li,v  that  might  exist  in  a 
narrow  region  of  slit  widths  and  system  interaction  parameters, 
satisfying  a  rather  subtle  balance  between  system  interactions  and 
block-wall  adsorption  energy.  The  weakened  influence  of  the 
nanopatterns  on  the  self-assembly  is  seen  mainly  for  W  >  16  and 
AflBM,2/A«AB  =  0-5,  where  the  systems  exhibit  chevron-like 
lamellar  phases  (Figure  13e).  In  contrast  to  the  A(7BM,2/AaAB  = 
0.5  case,  systems  with  wide  slit  widths  and  AaBMS/AaAB  =  1 
self-assemble  into  L— Li,2  +  iL — l^i  2,  which  nearly  follows  the 
superposition  behavior  of  the  parallel  alignment.  Additional 


Figure  12.  Self-assembly  of  the  A5B5  systems  confined  in  nano- 
patterned  planar  slits,  where  strips  on  opposing  surfaces  are 
vertically  staggered  with  size  Z,NP  =  12.5.  W is  the  slit  width,  and 
Aubh’j/Aoab  is  a  measure  of  phobicity  of  the  B  blocks  toward  the  w’2 
walls.  Schematics  indicate  lamellar  alignment  in  the  planar  slits 
without  nanopatterns  (i.e.,  in  the  planar  slits  with  uy  and  w2  walls). 
W]  and  u'2  walls  are  represented  by  black  and  blue  lines,  respec¬ 
tively.  Crosses  indicate  cases  where  nanostructures  did  not  result 
from  the  superposition  of  the  lamellar  alignment  in  the  underlying 
planar  slits  without  nanopatterns.  Letters  a— e  refer  to  simulation 
configurations  for  some  select  nanostructures  formed  in  the  nano- 
patterned  planar  slits,  which  are  shown  in  Figure  13.  Other  legend 
definitions  are  Lx,  a  perpendicular  lamellar  phase,  and  L||,2>  a 
parallel  lamellar  phase  with  two  A  lamellae.  For  the  case  in  which 
different  nanostructures  were  observed  for  differently  sized  simu¬ 
lation  cells,  both  nanostructures  are  shown,  separated  by  a  comma. 

DPD  simulations  with  box  lengths  of  Lx  =  25  and  Ly  =  20 
lead  to  a  different  nanostructure  only  for  the  case  of  [  W,  A a^wJ 
A<3Ab]  =  [8,  1]  (i.e.,  at  the  microphase-separation  boundary 
between  Lx  and  L— Li.i  +  iL— Li,i)>  where  the  system  formed 
Li, 2  rather  than  Lx.  (Both  nanostructures  are  indicated  in 
Figure  12,  separated  by  a  comma). 

The  schematic  phase  diagram  and  some  select  nanostructures 
for  narrow  (LNP  =  6.25),  vertically  staggered  nanopatterns 
together  with  a  detailed  discussion  are  left  for  the  Supporting 
Information.  Furthermore,  examples  of  conformational  behavior 
in  terms  of  the  nematic  order  parameter  for  the  A5B5  copolymers 
confined  in  the  vertically  staggered  nanopattern  slits  are  also 
provided  in  the  Supporting  Information. 

3.2.3.  Planar  Slits  with  Vertically  Partially  Staggered 
Nanopatterns.  For  vertically  partially  staggered  nanopatterns, 
the  self-assembly  behavior  does  not  follow  the  superposition 
description  for  either  the  narrow  or  the  wide  nanopatterns  except 
in  the  cases  of  very  narrow  slit  widths,  where  the  systems  form,  in 
addition  to  the  usual  perpendicular  lamellar  phase,  a  novel 
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Figure  13.  Examples  of  simulation  configurations  corresponding 
to  Figure  12.  (a)  L  +  iL,  (b)  L— Lil.i  +  iL — Ln  i ,  (c)  L— Lii,2  + 
iL — Lii,2,  (d)  a  U-shape  lamellar  phase,  and  (e)  a  chevron-like 
lamellar  phase.  Green  spheres  represent  A  beads,  and  gray  and 
pink  spheres  represent  w\- and  u’2-wall  beads,  respectively.  B  beads 
are  not  shown  for  visual  clarity. 

nanostructure  given  by  a  superposition  of  the  adsorbed  layer,  the 
inverted  adsorbed  layer,  and  the  perpendicular  lamellar  phases. 
For  wide  slits,  the  systems  typically  form  tilt-like  lamellar  phases. 
The  schematic  phase  diagram  and  some  select  nanostructures  for 
vertically  partially  staggered  nanopatterns  together  with  a  more 
detailed  discussion  are  left  for  the  Supporting  Information. 

3.3.  Connection  to  Real  Systems.  Thus  far  in  this  work,  all 
results  are  presented  in  reduced  units  with  no  connection  to  a  real 
system.  However,  our  DPD  model  of  a  symmetric  diblock 
copolymer  can  be  mapped  onto  actual  symmetric  diblock  copo¬ 
lymers  by  establishing  the  physical  values  of  the  cutoff  radius,  rc , 
the  mass  of  a  DPD  bead,  ;nDPD,  and  the  energy  term  kBT  by  using, 
for  example,  a  top-down  coarse-grained  approach.30  Closely 
following  the  work  in  ref  30,  we  consider  low-molecular-weight 
symmetric  diblock  copolymer  poly(styrene-6-2-vinylpyridine) 
(PS-PVP)  with  an  average  molecular  weight  of  10400  g  moF1 
that  consists  of  approximately  50  PS  and  50  PVP  monomers. 
Because  our  DPD  diblock  copolymer  model  corresponds  to  a 
flexible  Gaussian  chain  composed  of  Kuhn  segments,31  the 
number  of  monomers  defined  by  a  Kuhn  segment,  C„,  can  be 
estimated  using  Bicerano’s  structural-relation  method.32  Follow¬ 
ing  this  procedure,  we  determine  Co  —  9.9  for  both  the  PS  and 
PVP  blocks.  Hence,  a  PS-PVP  chain  of  50  PS  and  50  PVP 
monomers  can  be  mapped  onto  a  DPD  chain  of  5  A  and  5  B 
beads,  A5B5.  Defining  mmon  and  Vmnn  as  the  average  mass  and 
average  volume  of  the  PS-PVP  monomers,  respectively,  the  mass 
and  volume  of  a  DPD  bead  can  be  determined  from  mDPD  « 
CooWmon  -  1030  g  moF1  and  FDPD  C„Vmoa  -  1570  A3, 


(30)  Maly,  M.;  Posocco,  P.;  Pricl,  S.;  Fermeglia,  M.  Ind.  Eng.  Chem.  Res.  2008, 
47,  5023-5038. 

(31)  Rubinstein,  M.;  Colby,  R.  H.  Polymer  Physics ;  Oxford  University  Press: 
London,  2006. 

(32)  Bicerano,  J.  Prediction  of  Polymer  Properties;  Marcel  Dekker:  New  York, 
2002. 


respectively.  The  value  of  rc  can  be  defined  as  the  side  of  a  cube 
containing  an  average  number  of  p  beads: 

;-c  =  (pFdpd)1^3  (18) 

Equation  18  gives  for  the  PS-PVP  system  rc  =  16.8  A.  Using  the 
established  value  of  rc,  we  can  convert  the  geometry  of  the  planar 
slits  and  simulation  quantities  such  as  Rc2,  Rs2,  IVLAM,  and  WL 
into  real  units.  For  example,  the  slit  widths  in  reduced  units 
between  2  and  20  correspond  to  real  values  between  33.6  and  336 
A.  Similarly,  the  size  of  the  nanopattemed  strips  in  reduced  units, 
6.25  and  12.5,  correspond  to  real  values  105  and  210  A,  respec¬ 
tively,  which  can  be  compared  with  the  real  value  of  the  perpen¬ 
dicular-phase  lamellae  width  of  54  A  (about  3.2  in  reduced  units). 

The  value  of  kBT can  be  obtained  by  expressing  in  terms  of 

the  solubility  parameters  of  the  A(PS)  and  B(PVP)  beads,  (5a  and 


za£  =^((5a-(5b)2  09) 

For  the  PS-PVP  system,  Bicerano’s  structural-relation  method 
gives  dA  =  19.80  (J  cm3)1/2  and  dB  =  22.12  (J  cm3)1/2.  Therefore, 
when  yAH  =  2  is  used  in  eq  19  we  obtain  T  =  300  K  and  ksT  =  4  x 
IGF21  J.  To  complete  the  mapping  of  our  DPD  confined  diblock 
copolymer  system  onto  a  real  system,  one  needs  to  relate  our 
measures  of  phobicity  between  the  diblock  copolymer  beads  and 
the  walls,  AaAwJAaAB  and  AaBwJ A«ab,  with  the  real  interactions. 
In  contrast  to  pure  diblock  copolymer  systems,  where  mapping 
results  front  the  correspondence  between  DPD  and  Flory— 
Huggins-type  polymer  models  and  solubility  parameters,  it  is 
not  a  straightforward  task  to  map  the  diblock— wall  interactions. 
A  possible  route  is  by  matching  the  DPD  values  of  the  surface 
tension  for  the  block  polymer— wall  interfaces  with  the  corres¬ 
ponding  experimental  values.34  However,  we  leave  such  work  to 
another  study. 


4.  Conclusions 

We  performed  a  systematic  study  of  model  symmetric  diblock 
copolymers  A5B5  confined  in  planar  slits  with  and  without 
nanopatterns  using  DPD  simulations.  The  DPD  A5B5  copolymer 
can  be  considered  to  be  a  generic  model  for  low-molecular-weight 
symmetric  diblock  copolymers  such  as  the  PS-PVP  copolymer 
with  an  average  molecular  weight  of  1 0  400  g  moF 1  at  300  K.  The 
nanopatterns  were  modeled  by  alternating  portions  of  the  slit 
surface  that  interacted  differently  with  the  diblock  copolymer 
beads.  We  considered  various  nanopattemed  slits  where  patterns 
on  opposing  surfaces  were  vertically  (a)  aligned,  (b)  staggered, 
and  (c)  partially  staggered  for  two  differently  sized  nanopatterns: 
105  and  210  A  (or  6.25  and  12.5  in  reduced  units).  We  varied  the 
slit  width  from  33.6  to  336  A  (or  2  to  20  in  reduced  units)  and  the 
diblock  copolymer  bead— wall  interactions  from  preferential 
affinity  to  strong  phobicity.  We  determined  the  confined  phase 
behavior  by  analyzing  equilibrium  configurations  in  terms  of 
bead-density  profiles  and  isosurfaces,  along  with  the  coinciding 
characterization  of  the  confomiational  behavior  of  the  confined 
copolymers.  We  reiterate  that  the  goal  of  this  work  was  not  to 
determine  the  precise  location  of  the  microphase  separation 
boundaries  but  rather  to  investigate  systematically  and  qualify 
the  competing  effects  of  physical  adsorption  and  confinement  on 


(33)  Prausnitz,  J.  M.;  Lichtenthaler,  R.  N.;  de  Azevedo,  E.  G.  Molecular 
Thermodynamics  of  Fluid-Phase  Equilibria,  3rd  ed.,  Prentice-Hall:  New  York,  1999. 

(34)  Maiti,  A.:  McGrother,  S.  J.  Chem.  Phys.  2004, 120,  1594-1601. 
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copolymer  self-assembly,  which  to  our  knowledge  is  the  first  such 
study. 

For  the  planar  slits  without  nanopatterns,  the  systems  self- 
assembled  into  perpendicular  and  parallel  lamellar  phases  with 
different  numbers  of  lamellae  depending  on  the  slit  widths  and  the 
diblock  copolymer  bead— wall  interactions.  When  the  planar  slits 
were  composed  of  different  types  of  opposing  walls,  the  systems 
formed,  in  addition  to  the  expected  perpendicular  and  parallel 
lamellar  phases,  adsorbed  layer  and  adsorbed  layer— parallel 
lamellar  phases,  and  a  mixed  lamellar  phase.  The  observed  phase 
behavior  in  the  planar  slits  without  nanopatterns  closely  corre¬ 
lated  with  the  computed  conformational  characteristics  of  the 
diblock  copolymers  and  with  the  widths  of  lamellae  and  adsorbed 
layers.  For  example,  in  the  case  where  the  B  beads  had  a 
preferential  affinity  toward  the  slit  walls  the  transformation 
between  the  parallel  lamellar  phases  with  different  numbers  of 
lamellae  was  accompanied  by  an  elongation  of  the  copolymers  in 
the  perpendicular  direction,  along  with  an  increase  in  orienta¬ 
tional  order  and  an  increase  in  the  widths  of  lamellae  and  B-bead 
adsorbed  layers. 

The  self-assembled  behavior  observed  in  the  planar  slits  with¬ 
out  nanopatterns  enabled  us  to  rationalize  the  self-assembly  of  the 
A5B5  copolymers  in  the  nanopatterned  planar  slits  because  the 
nanopatterned  planar  slits  can  be  virtually  divided  into  portions 
that  correspond  to  planar  slits  without  nanopatterns.  We  inves¬ 
tigated  the  applicability  of  superimposing  the  observed  nano¬ 
structures  in  the  planar  slits  without  nanopatterns  to  predict  the 
nanostructures  formed  when  nanopatterns  were  present.  For  the 
planar  slits  with  vertically  aligned  nanopatterns,  the  resulting 
nanostructures  followed  the  superposition  description  quite  well 
for  wide  nanopatterns  and  rather  poorly  for  narrow  nanopat¬ 
terns.  We  observed  two  types  of  novel  nanostructures  given  by  (i) 
the  superposition  of  parallel  and  perpendicular  lamellar  align¬ 
ment  with  Schrek’s  surface  twist  boundary  between  the  parallel 
and  perpendicular  domains  and  (ii)  the  superposition  of  parallel 
and  inverse  parallel  lamellar  phases.  Similar  to  vertically  aligned 
nanopatterns,  the  self-assembled  nanostructures  in  the  planar  slits 


with  vertically  staggered  nanopatterns  obeyed  the  superposition 
description  reasonably  well  for  wide  nanopatterns  and  poorly  for 
narrow  nanopatterns.  Here,  the  systems  exhibited  two  types  of 
novel  nanostructures  given  by  (i)  the  superposition  of  adsorbed 
layer  and  inverted  adsorbed  layer  phases  and  (ii)  the  super¬ 
position  of  adsorbed  layer-parallel  and  inverted  adsorbed  layer- 
parallel  lamellar  phases.  The  system  also  formed  a  U-shaped 
lamellar  phase  in  a  narrow  region  of  slit  widths  and  system 
interaction  parameters.  For  the  planar  slits  with  vertically  par¬ 
tially  staggered  nanopatterns,  the  self-assembled  nanostructures 
did  not  follow  the  superposition  description  for  either  narrow  or 
wide  nanopatterns  except  for  very  narrow  slit  widths.  For  these 
narrow  slit  widths,  the  systems  formed,  besides  the  usual  perpen¬ 
dicular  lamellar  phase,  a  novel  nanostructure  given  by  the  super¬ 
position  of  the  adsorbed  layer,  inverted  adsorbed  layer,  and 
perpendicular  lamellar  phases. 
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